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We present a consistent framework of coupled classical and quantum dynamics. Our result allows 
us to overcome severe limitations of previous phenomenological approaches, like evolutions that 
do not preserve the positivity of quantum states or that allow to activate quantum nonlocality 
for superluminal signaling. A 'hybrid' quantum-classical density is introduced and its evolution 
equation derived. The implications and applications of our result are numerous: it incorporates the 
back-reaction of quantum on classical variables, it resolves fundamental problems encountered in 
standard mean field theories, and remarkably, also the quantum measurement process, i.e. the most 
controversial example of quantum-classical interaction is consistently described within our approach, 
leading to a theory of dynamical collapse. 
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Opinions vary about the coexistence of and the inter- 
action between classical and quantum systems. In ortho- 
dox quantum theory, classical macrosystems and quan- 
tized microsystems coexist, their interaction is described 
asymmetrically. The influence of macrovariables upon 
microsystems is precisely taken into account as exter- 
nal forces. The back-reaction of quantized microsystems 
upon classical macrosystems is largely ignored, except for 
detector variables which are typically sensitive to certain 
microvariables. The theory of this specific back-reaction, 
called measurement theory m, predicts the statistics of 
the final states after the interaction. Yet, the interpre- 
tation of quantized dynamics is exclusively based on this 
non-dynamical model of back-reaction (cf. collapse of 
the quantum state), without it we could not test the va- 
lidity of quantized dynamics at all. Possibly, quantiza- 
tion extends to macrosystems, indeed the criteria of being 
macro- or microscopic are loosely if ever defined. Con- 
trary to quantized microvariables, quantized macrovari- 
ables may have significant back-reactions on generic clas- 
sical macrovariables as well. This becomes apparent in 
the widely used mean-field approximation [E| which, how- 
ever, has several fundamental drawbacks [|| as we shall 
recall in this Letter. The measurement theory also de- 
scribes back-reaction, but only for idealized detectors. 
In some attempts to define quantum gravity, matter and 
some fields are quantized, while other fields (gravity in 
particular) are treated as classical, requiring thus a def- 
inite hybrid - i.e. a coupled classical and quantum - 
dynamics. Thus, a general model of the back-reaction is 
desirable. Such a theory would describe the " collapse" of 
the wave function dynamically Q, would replace mean- 
field approximations in a systematic way |&-0| , and could 
have deep implications for quantum cosmology. The 
first conceptual attempts |H-p| were followed by ups and 



downs [g| , until severe limitations have been clarified [^ . 
In this Letter we use a straightforward transformation 
of the problem which automatically leads to a consistent 
model of hybrid dynamics. 

The difficulties to overcome in our Letter are unre- 
lated to the high complexity of the emblematic mean- 
field equation of quantum cosmology. Actually, they 
lie in foundational principles. To illustrate these dif- 
ficulties, we assume a quantum Pauli-spin a interact- 
ing with a classical harmonic oscillator of Hamiltonian 
Hc{x,p) = i(p^ -f x^). The spin interacts with the "mag- 
netic" field of the oscillator via the Hamiltonian 



Hi = K&sp. 



(1) 



Regarding H = He + Hj as the total Hamiltonian, it 
is natural to prescribe the Heisenberg equation of mo- 
tion dtd' — ilHja] — iKp[a3,a] to the spin. The clas- 
sical oscillator momentum satisfies the Hamilton equa- 
tion dtp = —dxH — —X but the coordinate cannot sat- 
isfy dtx — dpH = p + Kas since x, being a real number, 
should not evolve into a matrix. The obvious way out 
is to replace the operator a^ by its quantum expecta- 
tion value, i.e. to apply some mean field approximation. 
Yet, if taken literally, this implies that quantum expec- 
tations can be deduced with arbitrary precision from the 
measurement of the classical variables x and p. Hence, 
the message is that the classical oscillator should, in 
some way, inherit quantum fluctuations from the spin. It 
comes to one's mind that x and p should be random vari- 
ables, but not arbitrary ones. As we shall demonstrate, 
mathematical consistency imposes that the classical vari- 
ables X or p must never take sharp values: The consistent 
theory assumes an unremovable coarse-graining |lQ,ni| . 
The mathematical issue is the following. The phase 
space density pc{x,p) of a classical canonical system C 



satisfies the Liouville equation of motion: 
dtPc{x,p) = {Hc{x,p),pc{x,p)}i 



(2) 



where Hc{x,p) is the Hamihon function and {f,g}p — 
dxfdpg — dpfdxg stand for the Poisson bracket. On the 
other hand, the density operator pQ of a (canonical, or 
maybe discrete) quantum system evolves according to the 
von Neumann equation: 



dtPQ = -■i[HQ,pQ] 



(3) 



where Hq is the Hamilton operator. To introduce inter- 
action between Q and C, we assume the hybrid Hamilto- 
nian in the following form: 



H{x,p) ^Hq + Hc{x,p) + Hi{x,p) 



(4) 



where the interaction term (and thus the total Hamilto- 
nian, too) is a Hermitian operator for Q, depending on 
the phase space coordinates of C. 

The mean-field approach assumes sharp classical co- 
ordinates xt^pt at each time and the current quantum 
expectation value HMF{x,p]t) = tr[H{x,p)pQ{t)] of the 
Hamiltonian is regarded as the effective Hamilton func- 
tion for the classical subsystem C. The coupled evolution 
equations then take this form: 



dtPQ{t) = ~i[H{xt,pt),PQ{t)\ 



(5) 



dtxt = dpHMF{xt,Pt), dtPt = ~dxHMF{xt,Pt)- (6) 

This approach has well-known deficiencies. In particu- 
lar, it gives no account of indeterminacy of the classical 
states x,p inherited from the quantum uncertainties of 
Pq. There is thus an essential nonlinearity in the mean- 
field von-Neumann Eq. (pi) which leads to fundamental 
conflicts with principles of locality 111]. Furthermore, 
the mean-field effective Hamilton function Hmf{x^p) 
will never be the proper representative of the interaction 
when quantum uncertainties in pQ are large. 

A promising conceptual approach [g uses the hybrid 
density p{x,p) to represent the state of the compos- 
ite system C x Q. If the subsystems are uncorrelated, 
then the hybrid density simply factorizes as pe{x,p)pQ- 
In the general case, the hybrid density p{x,p) should 
be an (a;, p)— dependent nonnegative operator, satisfy- 
ing an overall normalization condition tr J p{x,p)dxdp = 
1. One interprets the marginal distribution pc{x,p) = 
trp{x,p) as the phase space density of the classical sub- 
system C while the density operator pQ = j p{x,p)dxdp 
represents the unconditional state of the quantum sub- 
system Q. Conditional quantum states are natural to 
introduce at fixed canonical coordinates (a;,p) of the clas- 
sical subsystem Q : 



Pxp = p{x,p)/pc{x,p). 



(7) 



Aleksandrov g proposed the following evolution equa- 
tion for the hybrid density: 

dtp{x,p) = -i[H{x,p),p{x,p)] + (8) 

+ HH{x,p),p{x,p)}p - ^{p{x,p),H{x,p)}p . 

If Hi{x,p) = 0, then this equation splits into the stan- 
dard Eqs. (|) and (|). 

Let us test Aleksandrov's equation on the spin- 
oscillator system (p. The generic form of the hybrid 
state is 



Pix,p) = ^[1 + s{x,p)ff]pc{x,p), \s\ < 1, 



(9) 



where s{x,p) is the spin- vector correlated with the oscil- 
lator's state. We read off the conditional quantum state 
(0) of the spin: p^p — ^[l + s{x,p)a]. For the hybrid 
state (p| ) and interaction Hamiltonian (|l|), Aleksandrov's 
Eq. (^reads: 

dtp = (xdp - pdx)p - iKp[^3, p] - ^^[(73, dxp\+ . (10) 

This equation easily violates the positivity condition 
\s\ < 1 on p. E.g., the initial normalized polarization 
vector 



1 



-{2x,2p,x^+p^-l), 



(11) 



leads to |s| > 1 for all x > under the evolution ([iq). 
So, the naive Eq. (^) is inconsistent since it does not 
guarantee the positivity of the hybrid density p{x,p) [^. 
One sees that the mathematical textures of the classical 
(C) and quantum (Q) systems, though well understood 
each separately, are not at all trivial to couple. 

A royal road offers nonetheless. Let us quantize canon- 
ically C as well. We do so temporarily and, at the end of 
the day, we regard it classical again. The hybrid Hamil- 
tonian (m transforms into the total Hamilton operator 
of the fully quantized system C (g) Q: 



■.H{x,p):^ Hq+ ■.Hc{x,p): + ■.Hiix.p): 



(12) 



where :. . .: stand for normal ordering in terms of the usual 
annihilation and creation operators (x ± ip)/v2, respec- 
tively. Let the equation of motion for the total system's 
density operator p be the standard von Neumann one: 



dtp= -i[.H{x,p):,p] 



(13) 



Our royal road is based on coherent states m^ . Coherent 
states \x,p) are eigenstates of the annihilation operator: 



{x + ip) \x,p) — {x + ip) \x,p) 



(14) 



Using Bargmann's convention |T3[, the coherent states 
satisfy the following differential relation: 



{dx 
(dx 



idp) \x,p) 
idp) \x,p) 



{x-ip) \x,p) 
. 



(15) 
(16) 



The latter relation expresses the fact that the bras |a;,p) 
are entire analytic functions of the complex canonical 
variable x+ip, a crucial fact as we shall see. The coherent 
states form an overcomplete basis with normalization 

/ = J |a:,p)(x,p|H^^PHg!±P!))d^dp. (17) 

It follows from Eqs.(|l|Q that 

:/(a;,p): \x,p) = (18) 

X + ip + dx — idp p ~ ix + dp + idx 



:/ 



: \x,p) 



for an arbitrary normal ordered function of the quantized 
variables on the l.h.s. On the r.h.s., the same symbols 
:. . .: mean that all derivations must be done first! 

We apply a projection to the density operator p of the 
fully quantized system C iSi Q and thus reintroduce the 
hybrid density: 



p{x,p) =tr[{\x,p){x,p\ (g) I)p 



exp(-i(a:;^+p^)) 
27r 



(19) 



Indeed, this can formally be considered the hybrid den- 
sity of the composite system C x Q as if C were unquan- 
tized (i.e. classical) again. This is what we are going to 
do. We can thus derive the closed equation of motion 
of the hybrid density ln^) from the exact von Neumann 
equation (O) . Using the basic relation (O) and the iden- 
tity (16|), we obtain the desired evolution equation: 



dtp{x,p) 



-i:H X 



(20) 



dx + idr, , 9„ - id, 



jx n^ I'^'p 



V' 



p ^^x 



p{x,p)+H.C. 



This hybrid dynamic equation is our proposal to cou- 
ple classical systems to quantum ones canonically. Note 
that the hybrid density p(x,p) incorporates our statisti- 
cal knowledge of Q's quantum state, of C's classical state, 
and of their correlations. If Hj{x,p) = 0, then by inte- 
grating both sides of (^) over x and p we get the stan- 
dard von Neumann Eq. (|), as it should be. The trace of 
both sides, however, does not lead to the standard clas- 
sical dynamics (^. Instead, we get (for obvious reasons) 
the evolution equation of a Husimi- function mM . To low- 
est order in the derivatives, however, this is the classical 
Liouville Eq. (||). Hence classical dynamics is recovered 
if both the Hamilton function and the state distribution 
change slowly with x and p (see e.g. |15| and references 
therein) . 

Consistency of the hybrid equation of motion (|2Cl| ) is, 
contrary to the case of the naive Eq. (||), assured by con- 
struction. It preserves the positivity of the hybrid density 
p{x,p) along with a certain analyticity property. In fact, 
the projection (|l9|) leads always to hybrid densities of the 
form 



p{x,p) 



exp(-i(a;^-Hp^)) ^ 
27r ^ 



ipn{x - ip)Lp'^^{x + ip) 



(21) 



where (pn {x — ip) are unnormalized non-orthogonal state 
vectors for Q, being complex entire functions of the com- 
binations X — ip oi C's classical variables. This positive 
form is then preserved by our hybrid equation (po|). The 
analyticity condition (E^) restricts the possible hybrid 
states: sharp values of x,p and, generally, characteristic 
phase-space dependences inside single Planck cells are 
excluded. 

In particular, the "pure state" form of (pT|), i.e. with 
a single ipcp'' dyadic term, is also preserved. The hybrid 
Schrodinger equation of the hybrid state vector ip{x — ip) 
follows from (|2(]| ) |1£]. For completeness, we mention that 
the polarization (11) can be reproduced by ip{x — ip) — 
[{x - ip)\+) + \-)]/VS where agji) = ±|±). 

A post-mean-field approximation, incorporating some 
of the fluctuations of system Q into its back-reaction on 
the system C, is worth to derive. Consider the Isi-order 
expansion of the Eq. ( po| ) in the derivatives of H{x,p), 

dtp(x,p) = -i[H{x,p),p{x,p)] + 
+ i{H(x,p),p{x,p)}p- i{p{x,p),H{x,p)}p 
-^,[dxH{x,p),dxp{x,p)] - i[dpH{x,p),dpp{x,p)] . (22) 

This equation has two additional terms with respect to 
the (mathematically inconsistent) Eq. (j^). These two ad- 
ditional terms reduce the domain of inconsistency. More 
importantly, it can be shown that the above equation is 
equivalent to the exact Eq. IGu) if C is harmonic and its 
coupling to Q is linear in x and p. In particular, Eq. ( |2^ ) 
gives a mathematically consistent theory of fully quan- 
tized atomic systems (Q) interacting with the fully devel- 
oped c/assica/ radiation field (C). Moreover, its physics is 
equivalent with the true fully quantized radiation theory 

ITTI. 



By taking the trace of Eq. (22) one can show that the 
evolution of the classical states is a flow: 



dtx = {dpH{x,p))xp, dtp = -{dxH{x,p))^ 



(23) 



where 



stands for the expectation values in the 



current conditional quantum state pxp- Obviously this 
flow resembles locally the naive mean-field Eq. (g). Here, 
however, the classical state is randomly distributed: it in- 
herits the quantum fiuctuations of Q. On the other hand, 
the deterministic evolution of the state's distribution is a 
remarkable fact with significant consequences being dis- 
cussed elsewhere. 

Dynamical collapse of the quantum state is encoded in 
our hybrid evolution equation (EQ). It is a most spec- 
tacular feature. The "presence" of standard collapse will 
be demonstrated on the spin-oscillator model. For the 
hybrid state (pi), our evolution equation (|2C|) reads: 



^K[o-3,9p/5], (24) 



dtp= (xdp -pdx)p- 



which differs from the naive Eq. (|10|) by the presence 
of the Ath term on the r.h.s. This term guarantees 
the positivity of the hybrid state p{x,p) for all times. 
The oscillator plays the role of the Stern-Gerlach appa- 
ratus detecting the quantized spin-component a^. The 
pointer variable of the apparatus is x, it is set initially 
to zero with precision A = 1. Accordingly, we choose 
p'^{x,p) = (27r)~^ exp(— i(a;^ + p'^)) for the oscillator's 
initial state. We shall switch on the interaction Hamil- 
tonian at i = for a short time r compared to the oscil- 
lator period, keeping the effective coupling as strong as 
g — KT ^ 1. Actually, we replace k by g5{t) ||l|. Let us 
assume that the spin's initial state is the superposition 
IV') = c+|-|-) -f c_|— ) and the initial hybrid state ^ is 
the uncorrelated \'tp){tp\p'(,{x,p). Immediately after the 
interaction the hybrid state p{x,p) becomes 



\c+f\+){+\pl{x + g,p) 



->(-|Pc(^-.9,P) 



c+c*_\+){-\exp{-i9^ - i9P)Pci^,P) + H.C. 



(25) 



The trace yields the pointer's state distribution: 

pc{x,p)^\c+\^p^cix-g,p) + \c^\^p''cix + g,p) . (26) 

Since g ^ 1, we shall ignore the overlap between the two 
terms on the r.h.s., so we can say that the pointer x has 
swung out to g± A or to — (.g± A) with the probabilities 
predicted by the standard measurement theory. Invoking 
Eq. (m), we can easily read out the conditional quantum 
state of the spin from Eq. (pa) [since the off-diagonal 
terms are damped by exp{—^g^), we ignore them]: 



P.p = \±){±\ 



±g, p«0. 



(27) 



This shows the standard collapse of the spin's quantum 
state: the quantum state is correlated with the classical 
pointer's position. 

Discussing this Letter's results we repeat that we are 
aware of the ambiguous contemporary views concerning 
the concept of genuine hybrid systems. Nevertheless, the 
old Copenhagen interpretation as well as recent quan- 
tum gravity and quantum cosmological models assume 
such hybrid systems. We made the necessary compro- 
mises to neutralize strict no-go theorems. Our equation 
( P0| ) is a first example of hybrid dynamics which is both 
mathematically consistent and physically relevant. The 
applications of our equation are numerous, for example 
as phenomenological models whenever the mean-field ap- 
proximation is poor. Moreover, we derive the post-mean- 
field equations which describe back-reaction of quantum 
fluctuations to first order. On the foundational level, we 
point out that our hybrid dynamics ( pO[ ) reproduces the 
ideal quantum measurement, including the collapse of the 
wave function and the motion of the classical pointer. Let 



us also stress the close connection to current phenomeno- 
logical theories of dynamic collapse which follow from our 
hybrid dynamics |l6[ . Our hybrid theory is likely to be an 
integrating concept for treating quantum measurement 
dynamically and to overcome the inconsistent mean-field 
method in quantum cosmology. 
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